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Summary. — The discovery of the GRB-supernova association poses the question 
on the nature of the inner engine as the outcome of Type Ib/c supernovae. These 
events are believed to represent core-collapse of massive stars, probably in low- 
period stellar binaries and similar but not identical to the Type II event SN1987A. 
The branching ratio of Type Ib/c supernovae into GRB-supernovae has the remark- 
ably small value of less than 0.5%. These observational constraints point towards 
a rapidly rotating black hole formed at low probability with low kick velocity. The 
putative black hole hereby remains centered, and matures into a high-mass object 
with large rotational energy in angular momentum. As the MeV-neutrino emis- 
sions from SN1987A demonstrate, the most powerful probe of the inner workings 
of core-collapse events are radiation channels to which the remnant envelope is op- 
tically thin. We here discuss the prospect of gravitational-wave emissions powered 
by a rapidly rotating central black hole which, in contrast to MeV-neutrinos, can 
be probed to distances as large as lOOMpc through upcoming gravitational-wave 
detectors LIGO and Virgo. We identify the GRB- emissions, commonly attributed 
to ultrarelativistic baryon-poor ejecta, with a new process of linear acceleration of 
charged particles along the axis of rotation of a black hole in response to spin-orbit 
coupling. We include some preliminary numerical simulations on internal shocks 
produced by intermittent ejecta. The results showing a radial splash, which points 
towards low- luminosity and lower-energy radiation at large angles, possibly related 
to X-ray flashes. 

PACS 98.70.Rz - 98.62. 



1. Introduction 

The Type Ib/c GRB-supernova association, beautifully demonstrated by GRB030329/ 
SN2002dh, establishes SNIb/c as the parent population of cosmological gamma-ray bursts. 
This, in turn, confirms their association to massive stars proposed by Woosley |5 411 191 13*4] . 
Massive stars have characteristically short lifetimes, whereby GRB-supernovae track the 
cosmological event rate. This is quantitatively confirmed by agreement in their true-to- 
observed event rates of 450-500 based on two independent analysis, one on geometrical 
beaming factors and the other on locking to the cosmic star- formation rate (141 E7j . 
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Furthermore, Type II and Type Ib/c supernovae are believed to represent the endpoint 
of massive stars ^3^] in binaries, such as in the Type Il/Ib event SN1993J :29 . This 
binary association suggests a hierarchy, wherein hydrogen-rich, envelope retaining SNII 
are associated with wide binaries, while hydrogen-poor, envelope stripped SNIb and 
SNIc are associated with increasingly compact binaries ■ By tidal coupling, the 

primary star in the latter will rotate at approximately the obital period at the moment 
of core-collapse. With an evolved core these events are therefore believed to produce 
a spinning black hole |SH EJ El ■ 

GRB-supernovae are rare, in view of the observed small branching ratio of less than 
1% off their parent SNIb/c. While all SNIb/c may be producing black holes, only some 
are associated with the production of ultrarelativistic outflows powering gamma-ray emis- 
sions [3HI CHH • Quite generally, black holes produced in core-collapse possess random kick 
velocities of a few hundred km s _1 [7|. Quite generally, therefore, black holes formed in 
core-collapse are parametrized by their mass M, angular momentum J and kick velocity 
K. By its kick velocity, a newly formed black hole typically escapes the central high- 
density regions of the progenitor star before core-collapse is completed. This suggests 
an association of GRB-supernovae with the small sample of black holes possessing low 
kick velocities 05]. I n these cases, the black hole remains centered, surges into a high- 
mass object with a rotational energy of up to 66% of the maximal attainable value, and 
subsequently spins up by accretion or spins down in a suspended accretion state ■ 

Spacetime around a rotating black hole is described by the Kerr metric |2D]. It de- 
scribes an exact solution of frame-dragging induced by the angular momentum of a black 
hole. Frame-dragging appears in nonzero angular velocities of zero angular momentum 
observers. Observational evidence for frame-dragging rotating black holes is successfully 
pursued through X-ray spectroscopy by ASCA and XMM in surrounding accretion disks 
[411 1181 IT2j . Kerr black holes further possess an anomalously large energy reservoir in an- 
gular momentum. Per unit mass, this energy reservoir exceeds that of a rapidly spinning 
neutron star by at least an order of magnitude. Finally, this energy reservoir is entirely 
baryon-free. 

Frame-dragging creates a large differential angular velocity between the black hole, 
when rapidly spinning, and surrounding matter. When the latter is unmagnetized, the 
black hole continues to accrete matter, enlarging it and spinning it up towards an extreme 
Kerr black hole along a Bardeen trajectory Alternatively, when the surrounding 
matter is uniformly magnetized, corresponding to two counter-oriented current rings |50j . 
a spin-connection becomes effective between the angular momentum of the black hole 
and the inner regions of the surrounding accretion disk which forcefully spins down the 
black hole Especially at high spin rates, this process is facilitated by an equilibrium 
magnetic moment of the black hole in its lowest energy state. This process is self- 
sustaining, as most of the black-hole luminosity is incident onto the inner face of a 
surrounding torus, created with strong differential shear by competing torques acting 
on its inner and outer face. This energy input may be the agent driving a dynamo 
which creates superstrong magnetic fields. In a suspended accretion state, the large 
energy reservoir in black hole spin energy is released with an efficiency of r\ = Qt/Qh 
determined by the ratio of the angular velocity fix of the surrounding torus to the 
angular velocity of the black hole. This energy release, on the order of 10 53 erg. 
This output is larger than that seen at present in the electromagnetic radiation of GRB- 
emissions and supernova ejecta (about 10 51 erg, corrected for beaming and asphericities) . 
However, this additional energy output can be safely accounted for by "unseen" emissions 
in gravitational radiation and MeV-neutrino emissions 1501 0Hj ■ 
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Frame-dragging also acts along the axis of rotation. Here, it creates an energetic 
interaction along open magnetic flux-tubes ( "ergo-tubes" ) subtended by the event horizon 
of the black hole. These interactions define a new mechanism of accelerating magnetized 
ejecta to ultrarelativistic velocities. This spin-orbit coupling is a classical analogue of 
an earlier quantum mechanical treatment |4fi| . This application radically differs from 
the common view, that black-hole energetic processes are strictly confined to the action 
of frame-dragging in the ergosphere. We shall illustrate the case of intermittent ejecta 
with a numerical simulation, whose shocked interactions serve as input to baryon-poor 
outflows powering gamma-ray emissions |38l . 

2. Theoretical background 

An extreme Kerr black hole reaches an angular velocity of its horizon given by = 
1/2M. This angular velocity far surpasses that of surrounding matter, when separated 
from the event horizon by a distance on the order of the linear size of the black hole. 
According to the first law of thermodynamics, the efficiency of energy transfer from the 
black hole to the environment reaches a fraction 



of its rotational energy E rot = 2Msin 2 (A/4), where sin A the specific angular momentum 
of the black hole per unit mass and fir denotes the angular velocity of the surrounding 
torus. The torus itself is suspended by competing torques acting on its two faces, one 
positive torque on its inner face and a negative torque on its outer face due to generic 
losses in energy and angular momentum in magnetic winds. The question arises: what 
are the radiation channels through which the torus processes this input? We shall discuss 
this in the next section. 

The Kerr metric shows black hole spin creating frame-dragging, described by particu- 
lar contributions to the Riemann tensor jTHj. In turn, the Riemann tensor couples to spin 
|.'S2I l'S.'SI 1571 130| . Rotating black holes hereby couple to the angular momentum (and spin) 
of nearby particles. We note that the angular momentum per unit mass represents a rate 
of change of surface area, while the Riemann tensor has dimension cm~ 2 in geometrical 
units. Hence, curvature-spin coupling produces a force, whereby test particles follow 
non-geodesic trajectories |37| . These gravitational spin-angular momentum interactions 
are commonly referred to as gravitomagnetic effects by analogy to magnetic moment- 
magnetic moment interactions, even though there is an interesting difference in sign. To 
quantify this in the application to outflows from black holes, we shall focus on spin-orbit 
interactions along the axis of rotation following |4fi| . Using dimensional analysis once 
more, the gravitational potential for spin aligned interactions should satisfy 

(2) E = uJ, 

where uj refers to the frame-dragging angular velocity produced by the massive body and 
J is the angular momentum of the particle. Spinning bodies hereby couple to spinning 
bodies 

We shall detail the radiative process of Q and @ in case of a particular topology 
of the magnetosphere, consisting of a torus magnetosphere supported by two counter- 
oriented current rings and, formed by a changes of topology from the outer layers of the 
inner torus magnetosphere, an open magnetic flux tube along the axis of rotation [SHEDI- 
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3. Radiation channels of a non-axisymmetric torus 

The surrounding matter is strongly sheared by competing torques, which creates a 
torus with a super-Keplerian inner face and a sub-Keplerian outer face. For sufficiently 
strong viscosity - regulated by magnetohydrodynamical stresses - between the inner and 
the outer face, a state of suspended accretion can arise, in which the energy and angular 
momentum input from the black hole is balanced by energy and angular momentum 
losses to infinity. This suspended accretion state is described by balance equations of 
energy and momentum flux, and in many ways is equivalent to pulsars when viewed 
in poloidal topology |50| . The kinetic energy of the torus intruduces a bound on the 
poloidal magnetic field-energy that it can support. We recently derived an estimate for 
the a maximal ratio of poloidal magnetic field energy-to-kinetic energy of the torus of 
about 1/12 G2]|. 

Rapidly rotating black holes dissipate most of their spin-energy in the event horizon, 
set essentially by the spin-rate of the black hole. The lifetime of rapid spin hereby satisfies 

m 

This estimate is consistent with durations of tens of seconds of long gamma-ray bursts 

EH 

The gravitational wave-emissions due to quadrupole emissions by a mass-inhomogeneity 
SM T are 

(4) L gw = y (cuM) 10/3 F(e) ~ B (M H /R) 5 (SM T /M H f, 

where cj ~ M 1 / 2 /i? 3 / 2 denotes the orbital frequency of the torus with major radius R, 
M = (SM T M H ) 3/5 /(SM T + M H ) 1/5 ^ M h (SMt/M h ) 5 / 3 denotes the chirp mass, and 
F(e) denotes a geometric factor representing the ellipticity e of the orbital motion. The 
linearized result (J2J has been amply confirmed by approximately one-solar luminosity in 
gravitational waves emitted by PSR1913+16 (with ellipticity e = 0.62 [HI EH)- Here, 
we apply |0J to a non-axisymmetric torus around a black hole, whose mass-quadrupole 
inhomogeneity SMt is determined self-consistently in a state of suspended accretion 
for the lifetime of rapid spin of the black hole. Quadrupole mass-moments will appear 
spontaneously due to non-axisymmetric waves whenever the torus is sufficiently slender. 

The competing torques on the inner and outer face promote azimuthal shear in the 
torus, leading to a super-Keplerian state of the inner face and a sub-Keplerian state of the 
outer face. The torus becomes geometrically thick and may hereby become sufficiently 
slender to allow instability of quadrupole wave-modes m = 2 for a minor-to-major radius 
less than 0.3260 [101 • This produces gravitational radiation at close to twice the angular 
frequency of the torus. 

The equations of suspended accretion |50| can be solved algebraically, giving solutions 
for the total radiation energies emitted by the torus, expressed as fractions of the rota- 
tional energy of the black hole. Solving for the simultaneous output in in gravitational 
radiation, MeV-neutrino emissions and magnetic winds, we find 

,5, E ^0. 2 M a (£)(i|), f„ 
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in units of M Q = 2 x 10 54 erg. This appreciation J5J of GRBs typically surpasses the 
true energy E 1 ~ 3 x 10 50 erg in gamma-rays |14j by several orders of magnitude with 
potentially some exceptions Sub-dominant emissions in magnetic winds power an 

accompanying supernova at a factor rj less than that in gravitational radiation, 

(6) E „ = 4*10» etg (iL) 2 (^-). 

The energy in torus winds (jSJ) further provide a powerful agent towards various pro- 
cesses: an accompanying supernovae [8], possibly radiatively-driven and radio loud by 
dissipation of their magnetic energy |48| . collimation of the enclosed baryon-poor outflows 
from the black hole as well as a source of neutron pick-up by the same |27| . 

The energy output in thermal and McV-ncutrino emissions is a factor S less than that 
in gravitational radiation, or 



(7) E„ = 2 x 10 53 erg 



( M H 
0.30/ \7M a 



At this dissipation rate, the torus develops a temperature of a few MeV and produces 
baryon-rich winds ISU] ■ 

4. Spin-orbit interactions E = OJ 

A particle in a periodic orbit, e.g., a charge particle in a Landau state confined to a 
surface of constant magnetic flux, carries orbital angular momentum. Geometrically, this 
motion represents a rate of change of surface area an associated spacelike two-surface, 
following the dimensional analysis mentioned in the previous section. This two-surface 
can be identified by considering the helical motion as seen in four-dimensional spacetime 
|52| . The result is an interaction with the Riemann tensor - the Papapetrou interaction 
between spin or angular momentum with curvature |32l I33| . 

The Kerr metric embodies an exact solution to spin-induced curvature, expressed by 
the appearance of frame-dragging in the Ricmann-curvature. Not surprisingly, the Pa- 
papetrou spin-curvature coupling produces spin-coupling to the frame-dragging angular 
velocity. In a neighborhood of the axis of rotational of the black hole, this gives rise to 
(J2J). To sec this, we may use the non-zero components of the Riemann tensor of the Kerr 
metric in Boyer-Lindquist coordinates relative to tetrad 1-forms given by Chandrasekhar 
|10|. Accordingly, the Riemann tensor creates the radial force 

(8) F 2 = Ji? 3 i2o = J AD = -d 2 ujJ. 

The assertion follows from 



(9) E=r 

J r 



Fods = luJ. 



The result JHl also follows from a completely independent derivation, by considering the 
difference in total energy between two particles in counter rotating orbits about the axis 
of rotation of the black hole. We insist that these two particles have angular momenta of 
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opposite sign and equal magnitude, J± — g < p c / > u t (fl± - co), where u b denotes their velocity 
four- vectors. Thus, we have 

(10) J± = g^u^uj 2 - (g tt + (w*)- 2 )/g 00 = ±J. 

This shows that u* is the same for each particle. The total energy of the particles is given 
by E± = (m') _1 + fi±J±j and hence one-half their difference 

(11) E= l -{E + -E_)=LuJ. 

A perfectly conducting blob of charged particles in a magnetic field in electrostatic 
equilibrium is characterized by rigid rotation |43j with angular velocity fit,. In the frame 
of zero angular momentum observers, the local charge-density is given by the Goldreich- 
Julian charge density [lfij . In the lowest energy state with vanishing canonical angular 
momentum, the angular momentum of the charge particles satisfies J = eA<p, where 
e denotes the unit of electric charge and A<p the 4>— component of the electromagnetic 
vector potential A a 

Consider blob of charged particles in an open magnetic flux tube about the axis of 
rotation of the black hole with magnetic flux 2irA ( j ) . The number N(s) of particles per 
unit height s of the blob, therefore, satisfies N(s) = (fib — uS)A$/e, where e denotes the 
elementary charge. A pair of blobs in both directions along the spin-axis of scale height 
h hereby receives an energy 

(12) Euob = ojJNh = io(n b - u)A%h = (1 x 10 47 erg) B 15 h 3 M H, 

where Km — h/M denotes the linear dimension of the blob, B15 — B/10 15 . Here, we use 
the normalized function H = Alu (fob — tv^ , where lu — lj/Qh an d fib — flb/f^H- Charged 
particles in superstrong magnetic fields are essentially massless. The ejection of a pair of 
blobs with energy (|12|) hereby takes place in a light-crossing time of about 0.3ms for a 
stellar mass black hole, such as a seven solar mass black hole of linear dimension 10 7 cm. 
This corresponds to an instantaneous luminosity on the order of 3 x 10 50 erg/s. 

An intermittent source produces iterated emissions of blobs - "pancakes" in the ob- 
server's frame [SB] - which can dissipate their kinetic energy in collisions in the internal 
shock model of gamma-ray bursts |3S] OH] ■ Fig. (JIJ shows a relativistic simulation of 
such internal shock using JETLAB F90. 

This presentation in terms of blobs is complementary to discussions on continuous 
outflows |5S], wherein gamma-ray emissions are attributed to shocks due to by steepening 
by temporal fluctuations or intermittency at the source |5S], or late-time interaction with 
the environment. 

5. GRB-supernovae as long-duration burst sources for LIGO/Virgo 

The predicted gravitational radiation JSJ is in the range of sensitivity of the broad 
band detectors LIGO and Virgo 011130] shown in Fig. ©, as well as GEO [TT1E3 
and TAMA 0. Matched filtering gives a theoretical upper bound on the signal-to-noise 
ratio in the detection of the long bursts in gravitational radiation from GRB-supernovae, 
shown in Fig. (0 As a hydrodynamical source, the frequency will be unsteady at least 
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Fig. 1. - The head-on two magnetized ejecta ("kissing jets") in the comoving frame of the 
center of mass, representing the collision of a fast blob overtaking a slow blob at large distances 
from an intermittent source. The simulation by JETLAB F90 has incoming flows each with 
Lorentz factor 1.5, and densities a few times that of the unmagnetized environment. A high- 
energy density region forms about a comoving stagnation point with strong amplification of 
the transverse magnetic field, and a radial splash produced by a pressure-driven, subsonic and 
radially outgoing flow with accompanying cocoon. This transverse morphology points towards 
accompanying low-luminosity low-energy gamma-ray or X-ray emissions over large angles. 



on an intermittent timescales, so that a time-frequency trajectory method is probably 
more applicable in practice The signal-to- noise ratio will hereby be between first 

and second-order detection algorithms. 

As the outcome of Type Ib/c supernovae, GRB-SNe are an astrophysical source pop- 
ulation locked to the cosmic star-formation rate. Thus, their contribution to the stochas- 
tic background in gravitational radiation can be estimated semi- analytically, given the 
band-limited signals assuming B — A/ / f e of around 10%, where f e denotes the average 
gravitational-wave frequency in the comoving frame. Here, we can use the scaling rela- 
tions E gw = E M H /M , f e = f a M Q /M H where M = 7M , E = O.2O3M 77 O .i and 
fo = 455Hz?7o.i, assuming maximal spin-rates (E ro t — E ro t,max)- For non-extremal black 
holes, a commensurate reduction factor in energy output can be inserted. 

Assuming a uniform distribution of black hole masses, e.g., Mjj = (4 — 14) x Mq, 
uncorrelated to the relative angular velocities rj of the torus, summation gives an expected 
spectral energy-density @H| 

(13) < e' B (f) >= 1.08 x 10- 18 /b(x) erg cm" 3 Hz" 1 

where Jb(x) = /b {x)/xmxfB (•) is a normalized frequency distribution. The associated 
dimensionless amplitude \/<Sb(7) = \ // 2G/nc 3 f~ 1 F B 2 (/) , where F B = ce' B and G de- 
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Fig. 2. - The expected distribution of strain 
amplitudes of GRB-supernovae and the de- 
sign targets for the strain noise amplitude of 
the LIGO and Virgo detectors in the h(f)- 
diagram. The fiducial source distance is 
lOOMpc, corresponding event rate is about one 
per year. 
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Fig. 3. - The expected cosmological distri- 
bution of GRB-supernovae to the stochas- 
tic background in gravitational radiation, on 
the basis of their correlation to the cosmic 
star- formation rate N(z). (Reprinted from 
52 ©Cambridge University Press.) 



notes Newton's constant, satisfies 

(14) \/^(7) = 7.41 x W- 2 %ify\x) Hz" 1 / 2 

where fs(x) = fs(x)/maxfs(-), fs(%) = Jb{x)/x 2 , Likewise, we have for the spcctal 
closure density = fFB(f)/p c c 3 relative to the closure density p c = 3Hq/8itG 

(15) n B (f) = 1-60 x l(TVi/n(z), 

where fn(x) = fn(x)/maxfn(-), fn{ x ) = xfs{x) and Hq denotes the Hubble constant. 

The results are shown in Figs. (|2I3[) . The extremal value of Ob(/) is in the neigh- 
borhood of the location of maximal sensitivity of LIGO and Virgo. 
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